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Dukelsky, Martin-Delgado, Nishino and Sierra 0] (hereafter referred to as DMNS) inves- 
tigated the matrix product method (MPM) comparing it with the infinite-size density 
matrix renormalization group (DMRG) §]. For equivalent basis size, the MPM produces an 
improved variational energy over that produced by DMRG and, unlike DMRG, produces a 
translationally-invariant wavefunction. The DMRG results presented were significantly worse 
than the MPM, caused by a shallow bound state appearing at the join of the two DMRG 
blocks. They also suggested that the DMRG results can be improved by using an alternate 
superblock construction [B] • [B] for the last few steps of the calculation. 

In this comment, we show that the DMRG results presented by DMNS are in error and the 
artificial bound state produced by the standard superblock configuration is very small even 
for m = 2 states kept. In addition, we calculate explicitly the energy and wavefunction for 
the [B] • [B] superblock structure and verify that the energy coincides with that of the MPM, 
as conjectured in 0. 

The matrix product method allows full SU (2) symmetry to be utilized in a natural way, 
and is relatively easy to implement in a practical calculation. On the other hand, utilizing 
such non-Abelian symmetries in DMRG appears, at first sight, to be a more difficult problem 
and it was not until the invention of the interaction-round-a-face DMRG (IRF-DMRG) by 
Sierra and Nishino |Q that the first breakthrough occurred. Recently, it has been shown that 
non-Abelian symmetries can be integrated into the DMRG algorithm directly, without the 
need for a vertex-IRF transformation |^,^]. This results in several simplifications over the 
IRF-DMRG algorithm, while in principle the numerical results should be identical for the two 
methods. Table | shows the results of a re-examination of the spin-1 Heisenberg chain using 
the non-Abelian method, originally carried out in reference ||. The results of DMNS for the 
ground state energy determined by the MPM and the IRF-DMRG are listed as e MP DMNS and 
gDMRG dmns respectively. Current results using the superblock structure [£?]•[£>] and [£?]••[£>] 
are listed in columns e 1 SItc and e 2 sltcs respectively. The DMRG results of DMNS are claimed 
to use the [B] • »{B] structure, hence the energies in columns e DMRG DMNS and e 2 sltes should 
agree exactly. The discrepancy is, we believe, due to a problem with the implementation of 
the IRF-DMRG algorithm used by DMNS. While there is still a shallow bound state located 
in the center of the superblock, the effect on the ground state energy is several orders of 
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Table I - Energy density of the spin 1 Heisenberg chain as a function of number of states kept. 
e MPDMNS ande DMRGDMNS are f rom reference g . 



MP DMNS 



DMRG DMNS 



1 site 



„2 sites 



-1.333333 
-1.399659 
-1.401093 
-1.401380 
-1.401443 
-1.401474 



-1.333333 
-1.369077 
-1.392515 
-1.401380 
-1.401436 
-1.401468 



-1.3333333 
-1.3996590 
-1.4010933 
-1.4013806 
-1.4014447 
-1.4014757 



-1.3333333 
-1.3996237 
-1.4010886 
-1.4013798 
-1.4014430 
-1.4014756 



1.58x10^ 
4.06xl0~ 4 
5.39xl0 -5 
1.63xl0" 5 
7.77xl0 -6 
1.35xl0~ 6 



magnitude smaller than reckoned by DMNS. Figure [l] shows the bond energy (Si ■ S^+i) as a 
function of lattice position for to = 2 states kept. For the superblock structure [B] • [B], the 
bond energy is exactly translationally invariant as predicted by DMNS and is shown here with 
a solid line. In this case, the bond energy for different values of to (column e 1 sltc in table ||) 
agrees with that calculated by the matrix product method by DMNS (column e MP DMNS in 
table |), verifying the conjecture that the wavefunctions produced by these two algorithms 
coincide in the thermodynamic limit. 



Fig. 1 - Bond energy (Si ■ Si+i) as a function of lattice position for the DMRG wavefunction of the 
spin 1 Heisenberg chain with m = 2. 
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